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ABSTRACT 


Monitoring of environmental variables is a 
great challenge to derive reliable products that 
characterize the Earth's surface in an accurate 
manner. Land cover and its changes are 
valuable components in map making that 
focus on climate change studies and human 
impact on the environment. The North America 
Land Change Monitoring System (NALCMS) 
is an initiative which conducts land cover 
and change mapping for Mexico, United 
States and Canada. This trinational project 
1s coordinated by the Commission for 
Environmental Cooperation (CEC) which 
facilitates the work among countries. For 
Mexico, the change mask 2005-2010 has been 
recently completed and will be combined 
with similar national products of the other 
countries to detect changes across North 
America. This work describes how the mask 
for Mexico was generated by means of 
MODIS monthly composites from 2005 and 
2010, assessing biannual differences between 
spectral bands, NDVI and two  textural 
features. Landsat TM 3 images were also 
used to develop a training data set through 
visual change detection over central and 
northern portions of Mexico where most of 
the Landsat tiles were available. This visual- 
made training set was applied to calibrate the 
change detection algorithm. The paper shows 
the results for automated change detection and 
validation that focused on the correct 
delineation of detected change areas in order 
to assess final confidence. During validation, 
all polygons identified as change were 
assigned to six categories. In addition 
analysis of highlights case studies 1s presented. 


Keywords: MODIS, Land cover, Change 
detection, Change characterization, Spatial 
distribution. 


1. INPRODUCCION 


The coverage on the surface of the Earth and its 
changes have been widely studied since 1t has a 
main importance concerning to understanding 
the evolution of the societies, the development 
of models in order to study local and global 
scale phenomena, and monitoring how human 
activities impact on the planet. This sort of 
studies includes forest disturbance mapping 
(Masek et al. 2008, Pouliot et al. 2009), 
inventories (Feranec et al. 2010), resource 
management of protected areas (Fraser et al. 
2009) and planning (Manonmani and Suganya, 
2010). In this sense, the North American 
Land Change Monitoring System (NALCMS) 
provides information about the land cover at 
continental scale and is looking forward to 
developing annual land cover products and 
change detection maps between specific 
periods. NALCMS is a trinational initiative 
coordinated by the Commission for 
Environmental Cooperation (CEC) which 
facilitates projects between Mexico, United 
States and Canada. Since the land cover map of 
2005 was published in 2010, the project aims to 
develop a change mask between 2005 and 
2010 in order to update the initial 2005 map. 
As each country involved in the project 
develops its own change mask, for the 
Mexican portion the National Institute for 


Statistics and Geography (INEGI), the National 
Commission for Knowledge and Use of 
Biodiversity (CONABIO), and National 
Forestry Commission (CONAFOR) generated 
and assessed the national product which will 
be merged with US and Canada maps. The 
Mexican change map was obtained by an 
automated process based on monthly MODIS 
composites of bands 1 to 7 that were produced 
by the Canada Centre for Remote Sensing 
(CCRS) for the three countries (Latifovic et al. 
2012). Additional data such as the NDVI and 
texture filters were generated by CONABIO. 
The method applied for change detection 
regards  spectral differences between a 
comparable set of features from 2005 and 
2010, matching data from the same month. 
This 1s a common technique for change 
detection(Muchoney and Haak. 1994, Healey et 
al. 1994 among others). Only extreme 
differences between both years were considered 
to find the best change patterns. In order to 
choose the best lower and upper thresholds for 
these extreme differences, the values were 
divided into different combination of quantiles 
and amounts of features. These combinations 
were faced against a training data set developed 
with Landsat STM images over the central and 
northern portions of Mexico, where Landsat 
data are available for both years and cloud free 
data offer the best chance for visual land 
change identification. The comparison with 
training data produced different indices and 
the combination of features and quantiles 
which performed best was chosen as the final 
change mask for Mexico. 

Even though different studies have been 
developed in Mexico at local and regional scale 
(Berlanga-Robles and Ruiz-Luna, 2002, Mass, 
1999), just a few national efforts are available, 
(INEGI 2009, Mas et al. 2004, Couturier et al. 
2012). The main advantage of the NALCMS 
change mask 1s that 1t 1s generated by an 
automated algorithm which avoids human 
influence derived from visual analysis, besides 
this method allows detect changes over 
almost two millions of square kilometers in so 
much less time that required by human efforts. 
However, resulting change polygons must be 
reviewed in order to assess the accuracy of 
the final product and its confidence in different 
regions. This paper focuses on the outputs of 
validation that are presented in Colditz et al. 
(2012), and characterize the type of change as 


well depicts its distribution over Mexico. Three 
study cases are also analyzed due to the 
reliable impact they have on the land surface. 
Urban sprawl in Cancun, on the Yucatan 
peninsula in south- east (Veloz Avilés, 2011); 
and Monterrey in the northeast of Mexico can 
be related to fast economic growth in those 
areas. Water level changes at Alvaro Obregon 
reservolr in the northwest illustrate how the 
algorithm works and why 1ts outputs have to be 
validated. 


2. METHOD 


The process of change detection mask was 
based on  normalized monthly  composites 
ofMODIS images of band 1 to 7 (blue to short 
wave infrared) for the entire country, as well as 
the NDVI and textural features generated from 
normalized data. This process is based on the 
assumption that both years correspond 
temporally  (Colditz et al. 2012). The 
difference for each feature was computed 
between each composite from 2010 to its 
corresponding cone from 2005. In order to 
evaluate the best set of features and the most 
suitable thresholds of extreme values related to 
well-known change areas, many runs were tested 
and the outputs were compared with training 
data developed by visual analysis on Landsat5 
TM images. The most similar combination of 
inputs and difference values divided into 
quantiles were established as the final change 
mask. 


2.1Training data 

A set of LandsatS TM images of April and 
October from 2005 and 2010 was downloaded 
in order to evaluate changes during dry season 
and at the end of the rainy season, respectively. 
Only images for the central and northern 
portions of Mexico were available. The final 
composes of 74 pairs images of 133 required 
images to cover entire Mexico, Figure 1. 
Finally, areas evaluated as real change by an 
analyst were digitized 1f the polygon detected 
was persistent in both dry and rainy season. 
Thus the identification of land cover changes in 
the final reference map served to train the 
algorithm in order to find the best combination 
of features, lower and upper quantiles and 
number features aggregated to flag a pixel as 
change and get the best correspondence between 
both products. 


Landsat coverage 





2005 2010 





Water level 


Figure 1. a) LansatS TM images availability over the central portions of Mexico 1s shown in red, the black 
patch corresponds to area where no cloud free images were found. b) Polygon digitized as change 1s shown 
over a water reservolr, water lever difference 1s clearly visible 


2.2 Change detection 


Once the reference map was finalized, many 
attempts to find the best parameters ofbiannual 
difference were tested. If we consider that no 
change between both years means values close to 
zero, so the extreme values of the difference 
histogram could be considered as potential change. 
Besides aggregation of features was taking into 
account due to some land cover changes are not 
recorded by many spectral bands, thus the portion 
of features from the entire set which match the best 
with the reference map was established as 
frequency of change flagged in each pixel. 

In this sense, different runs considering lower and 
upper quantiles combinations demonstrated that 
the best result was obtained with 120 features, 
including 7 bands of each month, NDVI and 
two derived textural features (Sobel and 
Laplace) based on NDVI. The most appropriate 
lower and upper thresholds were 1 and 99% 
quantiles respectively, and the frequency of 
extreme values for those quantiles was 28. The 
output was compared to the reference map using 
the phi-coefficient. For the preliminary change 
mask, a phi of 0.2559 was reported. Further 
processing exclude areas well-known as no 
change due to  phenological stages in 
agricultural zones and smalls  patches of 
remaining clouds among others. 


The final processing step applied a minimum 
mapping unit (MMU) of 4 pixels in order to 
show only representative changes. The area 
detected as change 1s 780,462.5ha (0.39% of 
the total land surface) where well-known 
change areas were excluded, and 698,936.3ha 
(0.35%) with applying the MMU (Colditz et al. 
2012). The mask without MMU represents the 
maximum potential change over Mexico between 
2005 and 2010 in Figure 2. 


2.3 Validation 


Once the final change mask was delineated 
Figure 2, validation consisted about identify each 
single pixel inside polygons assessed by the 
algorithm. This process was conducted by an 
analyst work and classification as change or no 
change laid on his criterion. Visual analysis was 
made contrasting monthly MODIS composites 
generated with RGB combination of 721 spectral 
bands. As it was done with reference map, 
images from April at dry season and October at 
the end of the rainy season were employed for 
this analysis looking forward detecting 
persistence change in both cases for 2005 and 
2010. 
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Figure 2. Final change mask with maximal potential change, red polygons represents areas 
detected by the algorithm based on extreme biannual differences 


In order to re-digitize change polygons and 
assign change values to each pixel, an 11x11 
grid was generated surrounding analyzed 
polygons; this was helpful to identify 
missing elements omitted by the algorithm that 
needed to be connected to the identified change 
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patch. Other areas identified as change by the 
analyst but not detected by the algorithm were 
not evaluated as 1s shown in Figure 3. Values 
of O were assigned to no change and 1 to 
change. 


|] Change — by MODIS 


|] No change — by MODIS 


l  Change-— by analyst 
2 No change -— by analyst 


|] No change by MODIS but 
change by analyst — not 


assessed 





Figure 3. Criteria to evaluate change, red squares are area assessed by the algorithm; pixels with 1 
are change identified by the algorithm and the analyst, those marked with 2 in red squares are areas 
where the algorithm detected change but analyst disagrees. Pixels marked with 2 in white squares are 

no change surrounding areas. Blue pixels were not evaluated 


2.4 Change characterization 


During the validation process, change types 
were assigned with six main categories as 
these were possible to evaluate by visual 
analysis. The  polygons were  labeled as 
“water/land”, “vegetation loss”, “burnt area”, 
urban sprawl”, “no-change” and “other”, the 
most usual change classified as other was 
found on agricultural lands. It 1s important to 
mention that during the characterization 
process only the polygons assessed by the 
algorithm were evaluated. 

This information helps to understand how 
land cover evolves in different regions 
according to human activities as fast population 
growth in big cities and natural phenomena as 
wild fires or atypical precipitations as well. 


3. RESULTS AND DISCUSSION 


3.1 Accuracy assessment 

The overall agreement considering all 
change types as one single class except 
“no- change” shows 46.78%, and errors of 
commission and omission of 24.92% and 
44.61% respectively. 


The classes which perform best between 
reference and change mask excluding “no- 
change were “urban sprawl” with 
commission and omission errors of 11.56% 
and 21.81% respectively, and “deforestation” 
with 11.85% and 32.78% as well. This evident 
agreement in these classes relies on the 
assumption of the extreme spectral response 
of the land occupied by vegetation in contrast 
to Iimpervious surfaces. Remaining classes 
presented higher commission and omission 
errors as 1t is shown in Table 1. As 1t was 
expected, “no- change” class presented highest 
errors, mainly in commission (96.16%) 
although omission 1s also significant because, 
the analyst also assessed all areas surrounding 
change pixels. It has to be taken into account 
that our spatially-explicit assessment is much 
stricter than sampling- based attempts, because 
the area evaluated corresponds to each pixel of 
the change mask and connected pixels. 

In addition we have to consider the difference 
in resolution of both sources, 1.e. objects that 
were visually identified in Landsat cannot 
be detected with MODIS  coarse pixels. 
Besides, the human eye has a different 
perception when  detecting changes, for 
instance, by interpreting the surrounding area to 
find changes. 


Tablel. Percentage of commission and omission error for the six classes of change type, the 
corresponding area in relation to entire country is also included for each class 


Change type Commission (%) 


Water/Land 30.2637 


Area (Kmó) 


887.1 


Omission (%) Area (Km?) 


40.7954 1408.6 


Z Deforestation 11.8519 38 32.7784 137.8 
3 s 


Burnt area 0.9 0.625 69.9414 156.6 
4 Urban sprawl 11.5645 36.3 21.8195 17.5 
88.2 


SS No change 96.1653 


6 Other (Agriculture) 26.2391 


VTA 


193.6 


46.7683 


39.0635 348.9 





3.2 Change distribution 


Observed distribution over Mexico varies for 
each region; most patches are concentrated over 
central and northern portions of the country 
(Figure 4) and are mainly associated to 
“water/land” changes and “no change” class. In 
the south, no large patches were detected; 
probably due to change areas over this 
region are usually small and they could not 
be identified by the algorithm which is based 
on 6.25ha pixels. Besides, when we applied 
the MMU to the preliminary change mask no- 
significant areas are not mapped. 

In the Northwest, where main agricultural 
regions are located, there were 
polygons associated to water level changes 
in reservoirs which support irrigation in the 
coastal plain of the states of Sonora, Sinaloa 
and Nayarit. Also, important detected changes 
relate to variations in the water level 
observed in the Laguna “Ojo de Liebre” in 
Baja California Sur due salt concentration 
changes.  Additionally, near the city of 
Ensenada, Baja California, big patches of burnt 


so many 


were 





rt 


area were observed. Finally for this region, 
“urban sprawl” was evident in Tijuana over the 
border between Mexico and United States, where 
the city has been growing towards the South and 
East. 

The main concentration of no-change patches 

over northern-central part can be associated to 

soil humidity conditions in arid zones that 
represent different spectral response but which 
does not correspond to real land cover changes. 

Water level due to flooding 1s also evident 
in Figure 4. For instance, in Mexican 

southeastern plain areas in Tabasco state are 
annually flooding as consequence of summer 

extreme rains associated to hurricanes in 

many cases. Another example is observed in 

US- Mexican floodways; these divert water 
from Rio Grande to limit flooding in the 

Brownsville- Matamoros urban region. The 

big blue polygon over this area relate to an 

inundation that lasted for six months in 2010 

when category 2 hurricane Alex hit the coast 
of Tamaulipas June 394 2010, and a few days 

later followed tropical depression +2. Both 
events caused more than 500 mm of water 
landfall IBWC, 2010). 
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Figure 4. Change mask developed for entire Mexico, patches color indicate the type of change in 
accordance with 1x classes established to evaluate land cover transition 


Important “urban sprawl” was observed for 
many large cities, where migration from rural 
areas and fast economic growth in some 
cases offer conditions to urban expansion. 
Besides Tijuana in northwest, urban land 
cover changes are clearly observed around 
Monterrey, Nuevo León, mainly towards the 
North and East of the city. In Yucatan 
peninsula, changes were detected in urban 
areas of city of Mérida in the northern- 
central part where expansion does not seem 
to be directional showing patches all around 
the town. 

Probably the best example of land cover 
changes due to urban expansion and 
detected with MODIS is the city of Cancún 
at the northeastern coast of the Yucatan 
peninsula. Until the Mexican government 
decided to establish a planned touristic center 
in late 1970*s, Cancun was a small settlement 
of farmers of fishermen. In 2010 the 
population of the municipality of Benito 
Juárez increased to 661,176  inhabitants 
(INEGI, 2010). This caused an urban 
expansion rate of 616.4ha per year between 
2005amd 2010(Veloz Avilés, 2011). 


4. CONCLUSIONS 


This paper presented the process for 
automated change detection over Mexico 
based on the analysis of coarse resolution 
MODIS ¡images in order to generate a 
national change mask which can be merged 
with similar products of Canada and United 
States in the frame work of the NALCMS 
project. The final national products aim to 
present a 2010 updated version from 2005 
land cover map for entire North America 
region. 

The validation process of change mask, that 
instead of sampling evaluated each single 
pixel over the whole area, offered the chance 
to assign general change types according to 
observed features at 250m MODIS spatial 
resolution. Change types were classified into 6 
categories established on the capability to be 
visually identified. 

Final results allow determining 
commission and omission errors to entire 
change polygons set and being stratified by 
change type identified. 

The spatial distribution of change patches in 
relation to land cover transition, allowed 
identifying patterns related to different regions 
known natural processes of the landscape, and 
social elements as determinant factors. 
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ABSTRACT 


Incident solar light in terrestrial vegetation can 
be either reflected, absorbed or transmitted by a 
leaf. The absorbed light photons can be used 
for photosynthesis or can be dissipated into 
heat or emitted as fluorescence photons. In the 
latter case, photons are brought in an excited 
state of chlorophyll and emitted with a longer 
wavelength than light before excitation. The 
chlorophyll fluorescence is emitted, at leaf and 
canopy level, in the region 650-850 nm, with a 
peak located at 685-690 nm and another one at 
735- 745 nm. The emitted fluoresced radiance 
1s always super-imposed on the background 
reflected signal, but with a much lower 
intensity than the reflected radiance from the 
leaf, therefore 1t has been necessary to develop 
techniques that allow separating both signals 
from  remotely  sensed  radiance. This 
fluorescence signal is instantly related to the 
functional status of the  photosynthetic 
apparatus. Monitoring of this excess energy 
dissipation pathway provides a tool for 
detecting early stress symptoms of vegetation 
- that is, before damage 1s irreversible. This 
responsive feature of  fluorescence also 
supports practical applications such as 
agronomic crop improvement and forest 
vegetation assessment in the face of climatic 
instability. 

The growing interest on the potential of 
measuring  Cchlorophyll  fluorescence from 
remote sensing on a global scale recently led to 


the initiation of the Fluorescence Explorer 
(FLEX) mission. FLEX, currently under 
development phase A/B1 of the 8th European 
Space Agency (ESA) Earth Explorer Program, 
aims to be the first space mission primary 
dedicated to the estimation of chlorophyll 
fluorescence at a global scale. Our group is 
involved in various FLEX-related research 
activities and a brief overview of their latest 
status will be presented. 


Keywords: FLEX, Fluorescence, 
Photosynthesis, Chlorophyll 
1 FLUORESCENCE IN 


PHOTOSYNTHESIS 


Photosynthesis 1s the plant process that 
synthesizes organic compounds from 
inorganic substances in the presence of sunlight. 
The most important step, chemically, in 
photosynthesis is the conversion of carbon 
dioxide and water into carbohydrates and 
oxygen. The first photosynthesis part is an 
energy conservation phase that includes the 
formation of an electron transfer redox chain 
and the  consequent generation of a 
protons concentration  gradient. Thylakoids 
(sheet structures inside chloroplasts) contain 
protein complexes that perform the management 
of electrons and protons necessary for energy 
conversion. The protein complexes in thylakoid 
included are four; two of them are 


photosystems, namely, systems capable of 
absorbing and transforming energy photon in 
redox electrochemical energy. They are called 
photosystem I (PSD and photosystem II (PSID. 


chlorophyll-a molecule is not stable in that 
excited state and must return to its fundamental 
energy state, and therefore, must transfer the 
excitation energy. There are three ways in 


Each one has a different spectral absorption which the electron recovers 1ts fundamental 
shape of chlorophyll- a in the reaction centre. state: 1. Photochemical processes, 
The chlorophyll of photosystem I, called P700, photosynthetic activity, included  CO2 


absorbs light at infrared wavelengths (735-745 
nm) and the PSI chlorophyll called P680, 
absorbs at red wavelengths (650- 850 nm) 
(Figure la). The third complex is a connector 
between the two systems, called bgf cytochrome 
complex. Finally there is a fourth protein 
complex not involved in the transfer of 
electrons, but using the proton  gradient 
generated the other complex through the 
thylakoid to synthesize. When the chlorophyll-a 
molecule present in PSI and PSI is 1rradiated 
with electromagnetic radiation, 1s able to absorb 
a certain quantity of the same energy 
corresponding to a photon with the same energy 
quantity that is necessary to excite an electron 
from the molecule of chlorophyll-a and cause 
the jump of the electron to a higher 
energy orbital, leaving the chlorophyll-a 
molecule in an excited state (Figure 1b). The 
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assimilation, consuming the ATP and NADPH 
generated in the initial phase of photosynthesis. 
2. No photochemical processes, mainly the 
thermal dissipation (no radiant) of the excitation 
energy. 3. Light re-emission (fluorescence (E). 
As seen, the F signal is related to 
photosynthesis, allowing a direct estimate of the 
actual photosynthetic activity of the vegetation. 
In fact, F emission occurs in competence with 
photosynthesis and heat dissipation (Figure 1c). 
Under stress conditions, plants show a decrease 
of photosynthetic activity due to a change in 
the distribution of the energy absorbed in the 
three scope processes. However, F measurement 
1s difficult because 1t represents a small amount 
of the total radiance measured by a sensor to be 
masked by the light reflected by the plant itself 
(Figure 2) (Gomez-Chova et al., 2006). 





Figure 1. (a) The fate of light within leaves (ESA, 2008). (b) The different processes of light 
energy dissipation in photosynthesis (Cendrero, 2010). (c) Scheme of light-energy capture by chorophyll 
molecules and the different dissipation mechanisms used to return to the ground-state energy-level 
(ESA, 2008) 


The PSI and PSI are responsible for the F 
emission and as each of them has one shape of 
chlorophyll-a in the reaction centre 
spectrally different, therefore, the light re- 
emission 1s spectrally different. P700 absorbs 
light at infrared wavelengths and therefore is 


responsible for the F emission peak in the 
infrared (F740) and P680 absorbs at red 
wavelengths and therefore 1s responsible for the 
fluorescence emission in the red peak (1690) 
(Figure 2a). 
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Figure 2. Solar-induced Chl fluorescence and reflectance spectra of a sugar beet leaf. (a): F spectrum in 
the region between 650 and 800 nm. (b): The actual reflectance (black curve) was computed as the ratio 
of the reflected and incident fluxes; the apparent reflectance (red curve) was computed as the ratio of the 
total upwelling flux (reflection plus F) and the incident flux. The difference between the two quantities 1s 
highlighted in the grey boxes, showing the responses surrounding the F emission maxima (vertical dotted 
lines) (Meron1 et al., 2009) 


Although fluorescence, photosynthesis and non- 
photochemical protection are closely connected, 
the translation of  fluorescence data to 
photosynthesis 1s not trivial. In fact, while under 
low light unstressed conditions (no  non- 
photochemical protection mechanisms 
activated) a negative correlation exists between 
fluorescence and photochemistry, most studies 
have observed that in the presence of plant 
stress and high light conditions, fluorescence 
declines with photosynthesis as a result of 
protective mechanisms which take place in the 
leaf to prevent damage caused by harmful 
radicals formed in such stress conditions (Van 
der Tol et al., 2009). 


2. FLUORESCENCE AND REMOTE 
SENSING 


The study of solar induced fluorescence (SIF) 

(also known as “passive” fluorescence) 1s one 

important research objective for space-based 
platforms. In some ways 1t 1s a more mysterious 
phenomenon than actively-induced fluorescence 
because a given signal level cannot be 
unambiguously interpreted without 
consideration of aspects such as atmospheric 
and environmental information, or the pattern of 
SIF behavior over time (necessitating repeat 
measurements). SIF typically has been studied 
with respect to the changes in the ratio of the 


Red:NIR peaks that occur with physiological 
fluctuations. The ratio has been shown to be 
responsive to chlorophyll content, decreasing 
with increasing chlorophyll as a result of greater 
re- absorption of the red F in the leaf. SIF also 
responds to  photosynthetic  fluctuations 
independently of chlorophyll content 
(Lichtenthaler $ Rinderle, 1988). Comparisons 
between active and passive F measurements, 
photosynthesis yield, and photosynthesis rate 
have also been made (Carter et al., 1990; Moya 
et al., 2004; Meroni and Colombo, 2006). SIF 
has been successfully correlated to 
photosynthesis at the leaf scale, and in a 
preliminary way at the canopy scale. Simple 
ratios of the two peaks or other features such as 
the numerical value of the integral between 
peaks normalized over a parameter such as 
absorbed or incident PAR are examples of how 
SIF information may be considered for use in 
characterizing photosynthetic performance and 
resilience. 

The capacity to observe SIF presents a 
complete novel option for space-based remote 
sensing, by extracting information on the real- 
time workings of the photosynthetic machinery 
of plant canopies. Because the production of CF 
is related to the function of the two 
photosystems active during  photosynthetic 
initiation in green plants, 1t provides an actual 
elimpse into the dynamics of the photosynthetic 
process. As a photosynthesis indicator, F may 


be envisaged for use in: 

. tracking vegetation responses to stress, 
including  resilience and recovery from 
natural or human-induced perturbations 

> quantifying  photosynthetic efficiency, 
which can serve as a proxy for light use 
efficiency in gross primary productivity (GPP) 
and global carbon models 


directly helping to quantify GPP Plant F under 
solar illumination adds a weak signal to 
reflected solar radiation. If both fluorescence 
emission and surface reflectance are assumed to 
follow Lambert's law, the radiance upwelling 
from vegetation (L) at ground level is therefore 
composed of two coupled contributions, one 
reflected (r E/m) and the other emitted (F) 


> investigating implications for (Meron1 et al., 2009): 
photosynthetic rate, and possibly even more 

TÍAJE(A! 1 —J — A 
La) = HE + F(A) ¡Wm *sr "um P (1) 


where 12 is wavelength, r is reflectance (free of the emission component), and E is total irradiance incident 
on the target. Even though F contributes to the signal detected by a remote sensor, its magnitude 1s 
small (1-5% of the reflected radiation in the near infrared) thus making the decoupling of the two 
contributions difficult. Zarco-Tejada et al. (2000) first recognized the effect of F on reflectance and 
demonstrated that 1t is possible to detect the F signal using reflectance measurements. In fact, the 
reflectance factor usually computed by the RS community as the ratio between upwelling and incident 
fluxes 1s indeed polluted by the F contribution (Meroni et al., 2009): 
ML(A) MAA) 
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This quantity, referred here as apparent reflectance, r*, differs from the pure reflectance (r) as indicated 


by the right hand side of Eq. (2) and as shown in Figure 2b. 


The F signal is comparably stronger and can be 
detected passively in narrow dark lines of the 
solar and atmospheric spectrum in which 
Irradiance is strongly reduced (the so-called 
Fraunhofer lines). In the visible and near- 
infrared, the solar spectrum at ground level 
exhibits three main “Fraunhofer” features which 
have been exploited for F estimation: Ha due to 
hydrogen (H) absorption in the solar atmosphere 
(centered at 656.4 nm) and two telluric oxygen 
(02) absorption bands in the Earth atmosphere, 
O2-B (687.0 nm) and 02-A (760.4 nm). The 
incident at-surface irradiance at different 
spectral resolutions (0.005, 0.1 and 1.0 nm Full 
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Width at Half Maximum, FWHM) and at the 
spectral regions around these three absorption 
features 1s shown in Figure 3. Irradiance data 
were simulated by means of the Matrix 
Operator Model (MOMO) (Fell £ Fischer, 
2001). Among other inputs, the sun zenith angle 
was set to 30%, with a mid-latitude summer 
atmosphere, the surface at sea level, and the 
aerosol optical thickness at 550 nm equal to 0.2. 
Resampling to the different  spectral 
resolutions was performed assuming Gaussian- 
like spectral response functions with the given 
FWHM (Meron1 et al., 2009). 
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Figure 3. Incident irradiance (E) at ground level at the three absortion lines normally used for F 
retrieval and at three different spectral resolutions (Meroni et al., 2009) 
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Recently, alternative  retrieval  algorithms 
based on the spectral fitting of hyperspectral 
radiance observations, Spectral Fitting 
Methods (SFMs), have been proposed. FLD 
and SFMs were used to estimate F starting 
from Top Of Canopy (TOC) fluxes at very 
high spectral resolution (0.12 nm) and 
sampling interval (0.1 nm), exploiting the O2-B 
(687.0 nm) and 02-A (760.6 nm) atmospheric 
oxygen absorption hbands overlapping the 
fluorescence emissions at the red and far-red 
spectral window, resulting that violation of FLD 
assumptions results in a positive bias in F 
estimation at both oxygen absorption bands that 
cannot be avoided even at the high spectral 
resolution considered (Meroni et al., 2010). 


3. PAST PROJECTS 


Projects carried out over the past decade have 
served as valuable proof-of-concept studies for 
remote measurement of SIE: 

> SIFLEX (Solar Induced Fluorescence 
Experiment, executed over Finland in 
2003) demonstrated that SIF was discernible 
even in boreal conditions using tower-based or 
airborne instrumentation. 

> SEN2FLEX  (Sentinel-2 and FLEX 
Experiment in 2005) used a new airborne 
instrument AirFLEX and the CASI (Compact 
Airborne Spectrographic Imager) to reveal that 
the fluorescence signal contains information not 
apparent in classical remote-sensing 
observation. 

. CEFLES 2 (CarboEurope, FLEX and 
Sentinel-2, carried out in Southern France in 
2007) showed that canopy fluorescence 1s 
indeed closely correlated with ecosystem carbon 
uptake and that fluorescence data improve 
diurnal model predictions of GPP. On the 
regional scale, the fluorescence signal could 
be correlated to regional  airborne 
measurements of carbon fluxes. 

> FluorMOD project (Miller et al. 2005), 
was a seminal work funded by ESA, which 
undertook to model leaf and  canopy 
fluorescence. An integrated model was 
developed,  comprising the sub-models 
FluorMODleaf and FluorSAIL. The focus was 
to create a model for unstressed vegetation, 
though 1t did  accommodate  diurnal 
downregulation. The decision to emphasize 
fluorescence under benign conditions permitted 
understanding what “normal? fluorescence 
behavior looks like, so that one application of 


the FluorMOD results could be to help identify 
the existence of stress based on the deviation of 
measured results acquired from space against an 
unstressed “baseline” or “control”. FluorMOD 
predicted that SIF might, indeed, be measurable 
from a space-based platform. Developments 
linked to the FluorMOD work include an 
improved leaf emission spectrum taking into 
account emission of the two photosystems 
involved in the initial states of 
photosynthesis (Pedrós et al., 2008), and a 
further elucidation of FluorMODleaf that 
incorporated  re-absorption of  emitted 
fluorescence (Pedrós et al., 2010). This model 
simulates apparent spectral fluorescence yield of 
a fresh leaf between 640 and 840 nm for any 
excitation light of the visible spectrum. 

> In other recent work (separate from 
FluorMOD), a computational procedure was 
developed to resolve the respective 
contributions of PSI and PSI to the leaf 
chlorophyll fluorescence emission spectra at 
room temperature (Palombi et al., 2011). The 
approach used Principal Component Analysis 
of fluorescence emission in the O-I kinetic 
induction phase. A single- band spectrum for 
fluorescence emission peaked at about 735 nm 
for PSL and a 2-band spectrum with maxima at 
685 and 740 nm for PSIl was identified. A 
linear combination of only these two spectral 
shapes yielded a good fit for any measured 
emission spectrum of the leaf under 
investigation (Ficus benjamina) and may be 
useful for discerning the fluorescence emission 
contributions of the  photosystems under 
different conditions. Thus, the dynamics of 
energy distribution between the two 
photosystems may be  «monitored at 
physiological (-room) temperatures rather than 
having to use -196%C to study PSI emission as 
has been done in the past. Separation of the PSI 
and PSI emission components can improve our 
understanding of the changes in fluorescence 
signal induced by environmental factors or 
stress conditions, and may be helpful for 
interpreting remotely captured fluorescence. 


4. FLEX MISSION 


FLEX will be the first space mission 
designed to optimize remote observation of 
the chlorophyll  fluorescence in  terrestrial 
vegetation. The ESA mission will undertake to 
provide global maps of vegetation-fluorescence 
emissions. These might then be converted into 
photosynthetic efficiency for use in global 
monitoring of health and stress resilience. The 


Photosynthesis Study will provide useful 
recommendations for generation of spatial- 
temporal maps of vegetation health and stress 
status in near real-time. 

FLEX will take advantage of additional 
measurements provided by the Ocean and 
Land Colour Instrument (OLCI) and the Sea and 
Land Surface Temperature radiometer (SLST) 
onboard Sentinel-3. The combined information 
ls Critical to resolve the  disagreements 
between different land models that attempt to 
simulate future atmospheric CO2 concentrations 


through interactive terrestrial carbon balances. 
The objective 1s a global coverage mission, 
monitoring vegetation photosynthesis along 
the seasonal cycles and the 
activation/deactivation of the photosynthetic 
mechanisms where photochemistry 1s adjusted 
to a lesser or greater level of activity. Global 
coverage should be achieved for all land 
surfaces between -56% and +75 latitude 
including major islands. 

FLEX will be a satellite flying in tandem with 
Sentinel-3. Observations of FLEX and S-3 will 
be acquired within 6 to 15 s from each other. 
Payload is the FLuOkRescence Imaging 
Spectrometer (FLORIS) measuring 
fluorescence in O2-A and O2-B bands plus 
reflectance over red-edge and Photochemical 
Reflectance Index (PRD  spectral  ranges. 
FLORIS operates in pushbroom mode with a 
swath of 100 to 150 km (contained within one 
of the S-3 OLCI camera swath) and a pixel size 
of 0.3 x 0.3 km?. Data in four spectral bands: 
0O2-A (750-777 nm); 02-B (660-714 nm); red 
edge (714-750 nm); PRI (520-660 nm) with 
sampling of 0,1 nm in oxygen bands and 0,3 to 
2 nm in red edge and PRI bands calibration 
relies on both cross-calibration with the 
equivalent bands of OLCI and dedicated 
devices. FLEX and S-3 images co-registered 
by on-ground processing through correlation 
algorithms. Vega launcher as baseline and 
mission duration of nominally 3,5 years (to 
cover at least 3 full vegetation growth cycles) 
with goal of 5 years. 

The optimal spatial resolution for observing 
land targets and monitor surface processes 1s 
difficult to define, since almost all spatial 
scales are important.  Land-atmosphere 
interactions can hardly be modeled in high 
detail for large areas and in most cases the 
surface—atmosphere coupling is represented at 
reduced spatial resolution. In the particular 
case of carbon exchanges between surface and 
atmosphere, the combination of requirements 
provides an optimum spatial resolution in the 
range of 100-300 m, which is detailed enough 


to identify vegetated patches while resolving 
interactions at scales comparable to typical 
heights of the atmospheric boundary layer. Here 
a resolution of 100—-300 m would be 
reasonable as well and in line with other data 
sources and modeling efforts. 


5. BIOHYPE PROJECT 


Biomonitoring of natural vegetation allows 
the acquisition of well-defined samples at 
anaffordable cost and allows determining air 
pollution at different time-scales. It reflects 
longer-term changes of environmental quality, 
because plant leaves accumulate pollution over 
months, or even years for evergreen species. In 
the BIOHYPE project, biomonitoring  1s 
organized at three tree structural levels: (1) the 
sub-leaf level, (11) the leaf level and (111) the 
canopy level. An airborne hyperspectral based 
biomonitoring is applied to cover the goals at 
canopy scale. The overall research goal is to 
develop, test and  validate a  passive 
biomonitoring methodology based on 
hyperspectral observations and based on a field 
experimental approach which lapses over 
several spatial scales as well as over several 
tree structural and temporal levels, 
depending on the parameter considered. One 
of main goals of the project is the investigation 
at sub-leaf, leaf and canopy level of the spatial 
distribution of chlorophyll fluorescence. These 
objectives will be fulfilled by two broad field 
and flight campaigns in two cities of two 
different climatic zones: Ghent (Belgium) and 
Valencia (Spain). The scientific outcome of the 
project will foster thorough knowledge of the 
correlation between tree characteristics at the 
sub-leaf, leaf and canopy levels and air 
pollution. Moreover, a more  profound 
knowledge can be acquired concerning the 
correlation between the different tree level 
structural characteristics, and quite specifically 
the relevance and  biomonitoring potential 
of airborne  hyperspectral, as well as 
chlorophyll  fluorescence data. Up and 
downward fluorescence yield are, to some 
extent, linearly related to each other, but 
depending on the leaf anatomy and chlorophyll 
content, the ratio or fitting can differ. It 1s 
known that chlorophyll is one of the drivers of 
the first fluorescence peak, hereby reflecting the 
activity of  photosystem II (Figure 4) 
(SR/00/131, 2012). Results also show that SIF 
acts as a more sensitive indicator of (traffic) 
pollution stress compared to more static 
biochemical variables such as  chlorophyll 
content (Van Wittenberghe et al., submitted). 
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Figure 4. The fluorescence ratio F678/F741 upward and downward versus the estimated 
chlorophyll content 


6. EXAMPLE OF ESTIMATION OF 


SOLAR-INDUCED VEGETATION 
FLUORESCENCE FROM SPACE 
MEASUREMENTS 

In this example, a methodology for the 


estimation of  solar-induced  chlorophyll 
fluorescence from spaceborne and airborne 
sensors is presented. The fluorescence signal 1s 
included in an atmospheric radiative transfer 
scheme so that chlorophyll fluorescence and 
surface reflectance are retrieved consistently 
from the measured at-sensor radiance. This 
methodology is tested on images acquired by 





the Medium Resolution Imaging Spectrometer 
(MERIS) on board the ENVlIronmental 
SATellite (ENVISAT) taking advantage of its 
good characterization of the O2-A absorption 
band. Validation of MERIS-derived 
fluorescence is carried out by applying the 
method to data acquired by the Compact 
Airborne Spectrographic Imager (CASI-1500) 
sensor concurrently to MERIS acquisitions 
(Figure 5). CASI-derived fluorescence is in turn 
compared with eground-based fluorescence 
measurements, a correlation coefficient R* of 
0.85 being obtained (Guanter et al., 2007). 


Figure 5. Fluorescence maps obtained from a MERIS (a) and CASI (b) images (Guanter et al., 


2007) 
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RESUMEN 


En este trabajo se compararon valores de 
NDVI de Landsat 5 TM y de MODIS en dos 
campañas agrícolas  contrastantes en lo 
referido a las precipitaciones y se determinó la 
correlación entre datos de ambos sensores. El 
área de estudio se ubica al Oeste de Entre 
Ríos donde se identificaron, en cada campaña, 
diez lotes de maíz, de soja, de trigo y de 
soja de segunda. Los datos de precipitaciones 
pertenecen a siete estaciones meteorológicas de 
la Dirección de Hidráulica de Entre Ríos. 
Se realizó un análisis estadístico de los valores 
de NDVI teniendo en cuenta los períodos 
críticos durante los cuales se determina el 
rendimiento. En 2008 las precipitaciones 
fueron 52% menores al promedio histórico, 
en cambio, en 2009 resultaron un 20% mayor a 
dicho promedio. ¡Los valores de  NDVI 
presentaron diferencias significativas entre 
campañas, y en algunos casos -maíz, soja de 
segunda- según la fuente del NDVI, Landsat o 
MODIS. El máximo NDVI de Landsat se 
halló en maíz (0,89) y el mínimo en trigo 
(0,47); en MODIS, el máximo NDVI se 
observó en soja y soja de segunda (0,84) y el 
mínimo en trigo (0,46). Los valores máximos 
corresponden a la campaña 2009/10 y los 
mínimos a la campaña 2008/09. Los coeficientes 
de determinación hallados entre datos de 
Landsat y de MODIS resultaron elevados, con 
un mínimo R? de 0,73 y un máximo de 0,92 en 
soja. Esto determina que ambos sensores 
fueron capaces de detectar las diferencias de 
vigor de los cultivos en ambas campañas. 


Palabras clave: NDVI, Landsat, MODIS, 
cultivos, sequía. 


ABSTRACT 


In this study we compare Landsat 5 TM and 
MODIS NDVI values in two agricultural years 
with contrasting precipitation .amounts, and 
correlation between both sensors data was 
determined. The study area 1s located in 
western Entre Ríos, and we worked with ten 
fields of each crop under study (corn, 
soybean, wheat and second-crop soybean) in 
both years. Precipitation data came from 
seven meteorological stations that belongs to 
the Dirección de Hidráulica de Entre Ríos. 
We performed a statistic analysis of NDVI 
values considering critic periods for yield 
determination. In 2008, precipitations were 
52% less than historical average, but in 2009 
precipitations were 20% higher. NDVI values 
were significantly different comparing 
between years, and in some cases -corn and 
second-crop  soybean-  comparing between 
sources of NDVI, Landsat or MODIS. The 
highest Landsat NDVI value was observed in 
corn (0.89) and the lowest one in wheat 
(0.47); the highest MODIS NDVI value was 
found in soybean and in second-crop soybean 
(0.84) and the lowest one in wheat (0.46). 
Maximum values belong to 2009/10 and 
minimum values belong to 2008/09. High 
determination coefficients were found between 
Landsat and MODIS data, with the lowest R“ of 
0.73 in corn and the highest R“of 0.92 in 
soybean. This determines that both sensors were 
able to detect crop vigor differences between 
years. 


Keywords: NDVI, Landsat, MODIS, crops, 
drought. 


INTRODUCCIÓN 


La diferencia normalizada de energía reflejada 
en las regiones del infrarrojo y del rojo del 
espectro electromagnético ha sido el centro de 
atención para el desarrollo de índices 
cuantitativos con el objetivo de estimar la 
condición de la vegetación a través de imágenes 
de sensores remotos (Thiam y Eastman, 2006). 
Cuando la vegetación comienza a sufrir un 
estrés, esta condición puede ser revelada por 
cálculos de NDVI, debido a que el color verde 
típico de una vegetación sana y vigorosa 
comienza a tomar tintes marrones lo que 
aumenta la reflectancia en esas longitudes de 
onda. La vegetación estresada presenta una 
mayor reflectividad en la región del visible 
(400-700 nm) en comparación con la 
vegetación sana. Á su vez, este comportamiento 
se Invierte en la región del irc (750-1200 nm) 
y esta diferencia se puede utilizar para 
caracterizar el estado de los cultivos (Kyllo, 
2003). 

La sequía es una adversidad recurrente en la 
agronomía, afectando las cosechas agrícolas y 
la producción ganadera en distinto grado de 
intensidad. Los medios tradicionales que se 
utilizan para identificar, caracterizar o evaluar 
la sequía, se fundamentan en la aplicación de 
información meteorológica, expresando el 
desequilibrio hídrico a partir de indicadores 
elaborados con variables atmosféricas u 
obtenidos por medio de procedimientos de 
balance hidrológico. 

La precisión de procedimientos convencionales 
utilizados para evaluar la sequía decrece a 
medida que la proyección se aparta del lugar 
de origen de la información primitiva, lo que 
para las variables relacionadas con la 
precipitación torna el resultado incierto O 
tentativo a escasos kilómetros de distancia, 
incluso en regiones llanas (Hubbard, 1994; 
Rotondo y Seiler, 2001). 

Las mediciones satelitales para el seguimiento 
de la vegetación presentan, en cambio, las 
siguientes ventajas: son realizadas por un único 
observador, ofrecen una visión sinóptica, en 
muchos casos su uso tiene un costo nulo y 
alcanza una visión repetitiva de 


aproximadamente toda la superficie de la Tierra 
(Johnson et al., 1993), lo que constituye un 
recurso más idóneo para el seguimiento de la 
vegetación que los datos meteorológicos. 

Los trabajos de Ji y Peters (2003) utilizando 
series de tiempo de NDVI y su asociación 
con el Índice Estandarizado de Precipitación 
(Mc Kee et al., 1995) para detectar las 
regiones afectadas por sequías, y el de Zang et 
al. (2005) empleando el producto MODIS NDVI 
y MODIS LAI (índice de área foliar) para 
elaborar el índice de impacto de la variabilidad 
climática (CVIL, por sus siglas en inglés), que 
cuantifica el porcentaje de la producción anual, 
tanto ganado como perdido debido a la 
variabilidad climática, dan cuenta de la 
importancia de la teledetección en el monitoreo 
de sequías y tratan de cuantificar las pérdidas 
de producción agrícola causadas por ellas. 
También se ha simulado la condición y el 
rendimiento de cultivos utilizando datos 
de  Landsat y MODIS (Doraiswamy et al., 
2004). 

A su vez, el análisis multitemporal del NDVI 
es de gran importancia para determinar adelantos 
o retrasos en el desarrollo y la senescencia foliar 
y también ha sido utilizado para estimar 
distintos parámetros de los cultivos y separar 
aquellos con patrones fenológicos distintos 
(Fischer, 1994). 

El objetivo de este trabajo es analizar y 
comparar datos provenientes de los satélites 
Landsat 5 TM y MODIS en dos campañas 
agrícolas contrastantes en lo referido a las 
precipitaciones y determinar la correlación 
entre datos de NDVI de ambos sensores. 


MATERIALES Y MÉTODOS 


El área de estudio se encuentra delimitada 
por el paralelo 31%4133 S en el Norte y el 
paralelo 32*1158” S en el Sur, y desde el 
meridiano 60%1358” W al Este hasta el 
meridiano 60%W4112” W al Oeste. Estas 
coordenadas corresponden a la región centro 
Oeste de la provincia de Entre Ríos, 
Argentina (Figura 1). 
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Figura 1. Localización del área de estudio 


Para caracterizar las precipitaciones de la 
provincia de Entre Ríos durante las campañas 
agrícolas evaluadas (2008/09 y 2009/10), se 
utilizó el registro de datos de siete estaciones 
meteorológicas automáticas (EMA) de la 
Dirección de Hidráulica de Entre Ríos: 
Gualeguaychú, Villaguay, Basavilbaso, La Paz, 
Paraná, Puerto Yeruá y Diamante 
(http://www .hidraulica.gov.ar/precipitacionesH 
Istorico.php). 


Se localizaron mediante GPS diez lotes de 
cada cultivo evaluado en ambas campañas (trigo, 
maíz, soja y soja de segunda) y en cada 
lote se  confeccionaron 3 áreas de 
muestreo, cada una de 6x6 píxeles en las 
imágenes Landsat 5 TM (equivalente a 3,24 ha), 
totalizando 30 muestras por cultivo. En estos 
mismos lugares se registró el valor de NDVI 
proveniente de un píxel en las imágenes 
MODIS (equivalente a 6,25 ha). En la figura 2 
se puede apreciar un sector del área de 
estudio. 





Figura 2. Zona del área de estudio en donde se pueden apreciar algunas áreas de muestreo en una imagen 
A) Landsat 5 TM y B) MODIS 


Se trabajó con imágenes 
pertenecientes a las órbitas 226 y 227, fila 
MODIS/ Terra 
MODI13Q1, que ofrece información de índices 


82, y con el 


Landsat 


producto 


de vegetación integrando 16 días de datos. 


Para los meses en que normalmente ocurren los 
períodos críticos para la determinación del 
de los cultivos 
procedió a realizar un análisis 


rendimiento 


Tabla 1. Fechas de las imágenes utilizadas para comparar valores de NDVI entre campañas 


(Tabla 1), 
estadístico 


5 TM 


las imágenes 


Se 


Imágenes Landsat 5 TM 


Cultivo 2008/09 
Trigo 20-sep-08 
Maíz 16-dic-08 
Soja 26-ene-09 

Soja de segunda 22-mar-09 


satelitales 


para detectar diferencias entre los valores de 
NDVI de cada campaña y para comparar 
entre satélites ya que en los últimos estados 
vegetativos hasta inicio del llenado de granos, 
describen mejor el 
rendimiento esperado en varios cultivos (Yang 
et al., 2001). 


Imágenes MODIS 


2009/10 2008/2009 2009/10 

14-sep-09 13-sep-08 14-sep-09 
03-dic-09 18-dic-08 19-dic-09 
29-ene-10 02-feb-09 02-feb-10 
25-mar-10 22-mar-09 22-mar-10 


Las imágenes Landsat 5 TM se corrigieron 
atmosféricamente según el modelo propuesto por 
Chávez (1996) y se georreferenciaron con un 
error cuadrático medio menor a 15 m, utilizando 
una proyección Gauss-Kriiger (faja 5), datum 
POSGAR 94. El producto MODI13Q1 de 
MODIS posee una proyección sinusoidal. Los 
datos MODIS fueron obtenidos en línea de 
NASA Land Processes Distributed Active 
Archive Center (LP DAAC), USGS/Earth 
Resources Observation and Science (EROS) 
Center, SI0UX Falls, South Dakota 
(http://lpdaac.usgs.gov/get_data). 


RESULTADOS 


La campaña 2008/09 se caracterizó por una fuerte 
sequía que afectó a la región pampeana 
observándose una importante disminución de las 
precipitaciones. Promediando los datos de las 


siete EMA se puede apreciar que durante 
2008, sólo el mes de julio presentó una leve 
anomalía positiva en las precipitaciones con 
respecto al promedio histórico. En este año, se 
produjeron precipitaciones un 52% menor a dicho 
promedio. 

En cambio, en la campaña 2009/10 se registró 
una mayor cantidad de precipitaciones, 
aproximadamente 20% más que el promedio 
histórico durante 2009 (Figura 3). Entre 
septiembre de 2008 y enero de 2009 se registró 
una disminución de las precipitaciones algo 
mayor al 40 %. comparado con el 
promedio histórico para la zona. En cambio, en 
la campaña 2009/10 se registraron 
precipitaciones de aproximadamente un 70 % 
mayor para el período septiembre-marzo. 
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Figura 3. Anomalía de las precipitaciones en las campañas agrícolas 2008/09 y 2009/10 
(Elaboración propia con datos de la Dirección de Hidráulica de Entre Ríos) 


Vaiman (2011) encontró que la evolución temporal 
del NDVI de Landsat podía diferenciar el vigor 
de la vegetación entre campañas con marcadas 
diferencias en la cantidad de precipitaciones 


(Figura 4). Gran parte de estas variaciones 
en los valores de NDVI pueden explicarse por las 
variaciones de las condiciones climáticas 
(Jinyoung et al., 2010). 
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Figura 4. Evolución del NDVI en ambas campañas. Referencias: A: maíz; B: soja; C: trigo/soja 


Los valores de NDVI presentaron diferencias 
significativas entre campañas para ambos 
sensores, excepto para Landsat en soja de 
segunda (Tabla 2). Los valores mayores se 
obtuvieron een la campaña más húmeda 
(2009/10) en todos los cultivos. Los valores 
hallados en maíz en dicha campaña son 
superiores a los citados por Marini (2008) y 
Dediós et al. (2001). En soja, los valores de 
NDVI son similares a los mencionados por 
Jackson et al. (2004), quienes citan valores 
máximos de 0,9 en cultivos de soja con 


un elevado índice de área foliar y, en cambio, 
mencionan valores bajos de NDVI para 
cultivos con escaso desarrollo foliar. Esta 
situación es esperable por los efectos de una 
menor turgencia a nivel celular lo que afecta la 
expansión de las hojas cuando las plantas se 
someten a estrés hídrico (Pandey et al., 1984; 
Cox y Joliff, 1986). En trigo, Epinat et al. 
(2001) citan valores de NDVI de 0,92 para trigo 
en un año húmedo. Este valor es similar al 
hallado en este trabajo en trigo durante la 
campaña húmeda. 


Tabla 2. Comparación de valores de NDVI entre campañas para imágenes Landsat y MODIS 


Maíz 

NDVI Landsat* Soja 
Trigo 

Soja de segunda 
Maíz 

NDVI MODIS* Soja 
Trigo 

Soja de segunda 


Campaña 2008/09 Campaña 2009/10 
0,70 b 0,89 a 
0,54 b 0,87 a 
0,47 b 0,88 a 
0,75 a 0,79 a 
0,62 b 0,79 a 
0,52 b 0,84 a 
0,46 b 0,83 a 
0,66 b 0,84 a 


* Letras distintas en sentido de las filas indican diferencias significativas (p < 0,05; n = 30) 


Por otro lado, al comparar entre sensores 
para una misma campaña (p < 0,05; n = 30), 
los valores de NDVI en maíz y en soja de 
segunda presentaron diferencias significativas, 
en trigo existieron diferencias significativas solo 
en la campaña 2009/10, y en soja no se 
encontraron diferencias significativas entre 
ambos sensores. En todos los casos excepto 
en soja de segunda en la campaña 2009/10, los 
valores de NDVI provenientes de imágenes 
Landsat resultaron superiores a los de MODIS. 


Los valores de NDVI de  Landsat también 
resultaron superiores a los del sensor AVHRR 
de NOAA (Soria-Ruiz y Granados- Ramírez, 
2005). 

Los coeficientes de determinación hallados 
en la correlación entre los valores de NDVI 
del satélite Landsat y del MODIS resultaron 
elevados (Figuras 3 y 6). 
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Figura 5. Correlación entre valores de NDVI de Landsat y de MODIS. Referencias: A: maíz, B: 
soja; 1: campaña 2008/09, 2: campaña 2009/10 
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Figura 6. Correlación entre valores de NDVI de Landsat y de MODIS. Referencias: C: trigo, D: soja de 
segunda; 1: campaña 2008/09, 2: campaña 2009/10 


CONCLUSIONES 


Ambos sensores fueron capaces de detectar la 
diferencia en el vigor de los cultivos entre 
campañas, lo cual se explica por el déficit 
hídrico que sufrieron los cultivos en la 
campaña 2008/09 y las abundantes 
precipitaciones en la campaña 2009/10. 

Las diferencias encontradas entre sensores se 
podrían explicar por la diferente resolución 
espectral de cada uno y porque el producto 
MOD13Q1 representa una integración de 
datos de una serie de tiempo, y no un dato 
puntual como en el caso de las imágenes 
Landsat. 

La elevada correlación encontrada entre los 
datos de los dos sensores indicaría que 
ambos son adecuados para este tipo de estudio, 
pudiendo utilizarse uno u otro en base al 
nivel de detalle espacial requerido. Además, la 
mayor resolución espacial del sensor TM del 
Landsat discrimina con mayor eficiencia una 
cobertura de otra, y esto reduce la 
probabilidad de que se contamine la 
información con coberturas fronterizas. 
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ABSTRACT 


Chagas” disease whose incidence  1S 
conditioned to the presence of infected 
triatomines by Trypanosoma cruzi 1s a very big 
public health problem in the Brazilian Amazon 
region. Information on the eco- epidemiologic 
dynamics of vectors of diseases 1s needed to 
construct scenery of climate changes that have 
been noticed worldwide and related to 
antropic/local changes in the environment. The 
main propose of this work 1s to demonstrate the 
relationship between environmental 
degradation processes with the incidence of 
Chagas' disease in the  municipalities of 
Abaetetuba and Barcarena, in the period 
2000 to 2011. An environmental 
characterization using classified images of the 
Prodes and TerraClass Project. Abaetetuba 
showed a higher rate of deforestation and also 
an increased number of cases of Chagas 
disease. The rate of population increase, which 
was observed in the municipalities studied, too, 
can be a factor in increased deforestation and 
growth of cases of Chagas disease. Data from 
Prodes and TerraClass Project demonstrated 
the existence of a process of deforestation 
in the study area. The use of Geographic 
Information Systems with data from the 
Information System for Notifiable Diseases 
proved to be a very efficient tool to 
characterize the processes of environmental 
degradation that has occurred in the 
municipalities of Abaetetuba and Barcarena 
during the study period from 2000 to 2011. 


Palavras-chave: geographic information 
system, schistosomiasis, Biomphalaria, sistema 
de informacóes geográficas, esquistossomose, 
Biomphalaria 


1. INTRODUCAO 


Doenca de Chagas, cuja incidéncia está 
condicionada á  presenca de  triatomíneos 
infectados por Trypanosoma  cruzi é um 
problema de saúde pública muito grande na 
regiáo amazónica brasileira. 

Os vetores da doencade Chagas (DC) sáo 
insetos de orden  Hemiptera, subfamilia 
Triatominae e conhecidos popularmente no 
Brasil como  barbeiro.  kEstes  insetos 
possuem habito alimentar hematófago (Lent « 
Wigodzinsky, 1979; Sandoval et al., 2000). 

A doenca de Chagas e considerada uma doencga 
tropical, endémica em certas regi0es do Brasil e 
em outros 20 países da America Latina 
(Clayton, 2010). A doenca de Chagas afeta 
cerca de 10 milh0es de pessoas em países 
latino-americanos endémicos. No Brasil, 
estima-se que cerca de um milháo e novecentas 
mil pessoas possuam o mal de Chagas (Rass1- 
Junior et al., 2010). 

Na década de 60, foi relatado o primeiro surto 
agudo da doenca de Chagas na Amazónia 
brasileira. Fo1 descrito no estado do Pará e as 
investigacOes náo encontraram triatomíneos na 
área, propondo ser outra forma a via de 
transmissáo. Segundo Valente et al. (2009) 
cerca de 600 novos casos agudos da DC foram 
registrados na Amazónia brasileira. 


Nos últimos anos, vários surtos agudos da DC 
na Amazónia tém sido descritos, muitos desses 
associados a ingestáo do vinho do acaí (Euterpe 
oleracea) (Valente et al., 1999, 2009; Coura Jr., 
2006; Nobrega et al., 2009, Pereira et al., 
2009). 

Em funcáo destes surtos e outros demais na 
Amazónia, o Programa Nacional de Controle 
de Doenca de Chagas, da Secretaria de 
Vigilancia em Saúde (PNCDCh/SVS), vem 
implantando desde novembro de 2005, um 
modelo de vigilancia epidemiológica específico 
para a Regiáo Amazónica (Brasil, 2006). 
Estudos demonstram que a taxa anual de 
desmatamento na Amazónia brasileira e cerca 
de dois milhóes de hectares. De 1988 até 2010 
aproximadamente 385.234km” foram 
desmatados na Amazónia legal. Neste mesmo 
período, o estado do Pará foi o segundo estado 
com o maior índice de desmatamento no Brasil 
com uma taxa média de 5.543km”, perdendo 
apenas para o estado de Mato Grosso com 
5.834km” (Brasil, Projeto Prodes, 2012). 
Segundo Romana (2007) a transmissáo de 
doencas por vetores está associada com as 
mudancas da cobertura vegetal, particularmente 
com o desmatamento. 

O atual cenário da Amazónia brasileira 
configura no intenso fluxo migratório de 
pessoas procedentes de áreas endémicas, no 
desmatamento descontrolado, na grande 
diversidade encontrada de T.  cruzi, de 
hospedeiros vertebrados e  espécies de 
triatomíneos (Aguilar et al., 2007). 

O uso de ferramentas de geoprocessamento 
vem sendo amplamente utilizado em estudos na 
área de saúde, especialmente em pesquisas 
relacionadas com doencas tropicais. Com o uso 
destas ferramentas é possível fazer uma 
análise da  distribuicáo geográfica das 
espécies, usando modelos de nichos 
ecológicos, assim como, avaliar o risco de 
infeccáo ás populacoes e o seu controle. 
Sistemas de Informacáo Geográfica (SIG) sáo 
úteis nos estudos de doencas. Em geral, o 


desmatamento de grandes áreas  rurals 
interferirem com o potencial vetor, infeccáo de 
T. cruzi e de reservatórios de mamíferos, e que 
sáo desconhecidos para a regiáo ao longo dos 
anos.  Informacóes sobre a dinámica 
ecoepidemiológica de vetores de doencas sáo 
necessário para construir cenários de mudancas 
relacionadas ás alteracóes antrópicas. Desta 
forma, é possível avaliar as  condicóes 
favoráveis para a dispersáo dos vetores, bem 
como o seu desenvolvimento, o  sucesso 
reprodutivo e a sobrevivéncia. 

A principal proposta deste  estudo  é 
demonstrar a relacáo entre os processos de 
degradacáo ambiental com a incidéncia da 
doenca de Chagas nos  municípios de 
Abaetetuba e Barcarena, no período de 2000 a 
2011. 


2. METODOLOGIA DE TRABALHO 


A área de estudo foi os municípios de 
Abaetetuba que localiza-se a uma latitude 
01%43'05”5 e a uma longitude 4852*57”0, 
tem uma populacáio de 147.846 habitantes em 
uma área de 1.617 km” e o município de 
Barcarena que localiza-se a uma latitude de 
em uma área de 1.315 km” (IBGE, 2012). 
Ambos os municípios estáo localizados no 
nordeste do Estado do Pará e possui um clima 
equatorial super úmido (Figura 1). 

Foram utilizadas imagens do satélite Landsat 
223/61 e 224/61 (orbita/ponto) de 2000 a 
2011, obtidas do Projeto Prodes 
(http://www.obt.inpe.br/prodes/index.html) e do 
Projeto Terra Class 
(http://www.inpe.br/cra/terraclass.php), para 
fazer O mapa da cobertura vegetal dos 
municípios estudados e comparar com os 
casos de Doenca de Chagas obtidos da base 
de dados dos casos confirmados notificados 
no Sistema de Informacáo de Agravos de 
Notificacgáo (SINAN). Foram obtidos 136 casos 
de pacientes positivos para a Doenca de Chagas. 
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Figura 1. Localizacáo espacial dos municípios de Abaetetuba e Barcarena 


Foi desenvolvido um banco de dados com 
dados ecoepidemiológicos usando Sistemas de 
Informacoes Geográficas. Após o processo do 
banco de dados foram realizadas consultas 
espaciais, ligando a classificacáo do Projeto 
Prodes/TerraClass com os dados de Doenca de 
Chagas nos municípios de Abaetetuba e 
Barcarena. 

O software ArcGis (versáo 9.3 — 


http://www.esri.com/) foi utilizado para as 
analises espaciais. 


3. RESULTADOS E DISCUSSAO 


A Tabela 1 apresenta a quantidade de casos 
de pacientes positivos para a Doenca de 
Chagas nos municípios de Abaetetuba e 
Barcarena, no período de 2000 a 2011. 


Tabela 1. Casos de pacientes positivos para a Doenca de Chagas no período de 2000 a 2011 


Município 
Abaetetuba 0 0 1 1 
Barcarena 0 0 0 0 


Observa-se na Tabela 1 que nos anos de 2000 e 
2001 náo houve casos de pacientes positivos 
para a Doenca de Chagas, segundo dados do 
SINAN. Por isso, foram retirados das análises 


2000 [001 [002 [2003 [004 [P005 P006 R2007 pP008 pP009 2010 P011 


0 13 1 10 55 7 9 
0 1 1 ' 4 0 


esses dois anos, assim como, náo foram 
utilizadas as imagens classificadas do Projeto 
Prodes. A imagem classificada do Projeto 
Prodes do ano de 2011 pode ser visualizada 
na Figura 2. 
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Figura 2. Imagens classificadas do projeto Prodes dos municípios de Abaetetuba e Barcarena nos anos 
de: (a) Classificacáo do Projeto Prodes e (b) Classificacáo do Projeto TerraClass 


Dados dos Projetos Prodes (Figura 2A) e 
TerraClass (Figura 2B)  demonstram a 
existéncia de um processo de desmatamento na 
regiáo de estudo. A Figura 2(A) mostra o mapa 
de desmatamento da área de estudo do ano de 
2011 com os incrementos de desmatamento dos 
anos de 2001 a 2011. Observa-se na Figura 
2(A) que existe um processo de degradacáo 
ambiental mais elevado no :município de 
Abaetetuba do que em Barcarena. 

A Figura 2(B) mostra a classificacáo do 
Projeto TerraClass do ano de 2008, onde 
observa-se que  Abaetetuba e  Barcarena 
apresentam uma grande área agrícola. Também, 
nota- se que Abaetetuba tem uma grande área 


de vegetacáo secundária, além de várias áreas 
de regeneracáo com pasto. 

Ao longo da história de  análise de 
desmatamento da Amazónia Legal pelo Projeto 
Prodes (Brasil, Projeto Prodes, 2012), a 
quantidade de desmatamento de Abaetetuba foi 
de 42,68%, o que representa quase o dobro do 
ocorrido no município de Barcarena que foi de 
24,11%, além do fato de Abaetetuba possuir 
uma populacáio de 141.100 habitantes, que 
representa cerca de 40% a mais do que 
Barcarena que tem 99.859 habitantes (IBGE 
2012), esse fato pode ter contribuído para o 
processo de transmissáo da DC. 


Tabela 2. Porcentagem de desmatamento obtido do Projeto Prodes nos municípios estudados no 
período de 2002 a 2011 


Município 


2002 [2003 [004 [2005 |(P006 [007 pP008 (2009 pP010 2011 





Abaetetuba 0,16 (0,15 (0,7 0,2 
Barcarena 0.56 (0,75 0,69 [0,29 


Por outro lado, na Tabela 2 observa-se que no 
período de estudo (2002 a 2011) ao contrário 
do processo de desmatamento citado no 
parágrafo acima, Barcarena sofreu um maior 
desmatamento do que Abaetetuba, cerca de 
3,6% e 2,3%, respectivamente, dados esses 
também obtidos do Projeto Prodes (Brasil, 


0,28 
0:32 


0,25 0,08 0,19 [0,24 [0,05 
0,25 0,09 [0,21 0,14 [0,08 


Projeto Prodes, 2012). Este fato pode estar 
relacionado com a criagáo do polo industrial de 
Barcarena e, também, segundo o IBGE (2012) 
a populacáo de Abaetetuba cresceu cerca de 
20% nos últimos 10 anos (censo 2000 e 2010) e 
Barcarena cresceu cerca de 57%. 


4. CONCLUSOES 


A utilizacgio dos Sistemas de  Informacáo 
Geográfica (SIG) com os dados obtidos do 
Sistema de  Informacáo de  Agravos de 
Notificacáao (SINAN)  mostrou ser uma 
ferramenta muito eficiente para caracterizar Os 
processos de degradacáo ambiental que ocorreu 
nos municípios de Abaetetuba eBarcarena, no 
período estudado de 2000 a 2011. Também, o 
SIG demonstrou ser um poderoso recurso para 
captar, armazenar e recuperar informacóes. 

A classificacáo utilizando imagens digitais 
possibilitou visualizar nas áreas de estudos as 
características  ambientais  propícias a 
propagacáo da doenca de Chagas. 

Os resultados obtidos do Projeto Prodes e do 
IBGE sugerem a possibilidade que a ocorréncia 
da DC tenha sido influenciada por alteracóOes 
induzidas pelo homem no mejo ambiente. 
Também, os resultados mostraram que o 
município de Abaetetuba apresentou maior 
desmatamento do que no  município de 
Barcarena e isso pode ser a causa do aumento 
de casos de doenca de Chagas apresentados em 
Abaetetuba. Dias (2011) estudou a transmissáo 
de doenca de Chagas nas localidades de Sáo 
Tomé, Araipá e Nova Estrela pertencentes ao 
município de Aveiro/PA e também mostrou a 
relacáo entre desmatamento e DC. 

A taxa de crescimento populacional, que foi 
observada nos municípios estudados, 
também, pode ser um fator para O aumento 
do desmatamento e do crescimento de casos 
de doenca de Chagas. 

Seráo obtidos os dados de DC do banco de 
dados do IEC para uma  comparacáo, 
validacáo e, também, para a obtencáo das 
coordenadas das localidades/enderecos dos 
pacientes positivos para a Doenca de Chagas. 
Essas coordenadas seráo  georreferenciadas 
com o uso de um aparelho GPS e depois 
analisadas com os dados das imagens do 
Projeto Prodes para uma melhor identificacáo 
do tipo de uso do solo associado com a 
incidéncia da doenca. 
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